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^^' Using repeated Laplace transform techniques, along with newly- developed accurate numerical 

C^ I inverse Laplace transform algorithms [ll, Ql, we transform the coupled, integral- differential NLO sin- 

glet DGLAP equations first into coupled differential equations, then into coupled algebraic equations, 
which we can solve iteratively. After Laplace inverting the algebraic solution analytically, we numer- 
ically invert the solutions of the decoupled differential equations. Finally, we arrive at the decoupled 
NLO evolved sokitions 






Fs{x,Q^) = Ts{Fso{x),Go{x)), 
G{x,Q^) = g{Fso{x),Goix)), 

where J^s and Q are known functions- determined using the DGLAP splitting functions up to NLO 

in the strong coupling constant aa{Q^)- The functions Fso{x) = Fs{x,Qq) and Go{x) = G(x,Qq) 

are the starting functions for the evolution at Q^ = Qq. This approach furnishes us with a new tool 

for readily obtaining, independently, the effects of the starting functions on either the evolved gluon 

or singlet structure functions, as a function of both Q^ and Qq. It is not necessary to evolve coupled 

integral-difi'erential equations numerically on a two-dimensional grid, as is currently done. The 

^ ' same approach can be used for NLO non-singlet distributions where it is simpler, only requiring one 

>0 , Laplace transform. We make successful NLO numerical comparisons to two non-singlet distributions, 

"<) • using NLO quark distributions published by the MSTW collaboration [31, over a large range of x 

**! ' and Q^. Our method is readily generalized to higher orders in the strong coupling constant as{Q^ 

O ' I- INTRODUCTION 

o' 

In order to interpret the experimental results at the Large Hadron Collider — in the search for new physics — 
accurate knowledge of gluon distribution functions at small Bjorken x and large virtuality Q^ plays a vital role in 
ITJi I estimating QCD backgrounds and in calculating gluon-initiatcd processes. The gluon and quark distribution functions 
'V^ ' have traditionally been determined simultaneously by fitting experimental data on neutral- and chargcd-currcnt deep 
d ' inelastic scattering processes and some jet data over a large domain of values of x and Q^. The distributions at small 
X and large Q^ are determined mainly by the proton structure function i^^^(a;, Q^) measured in deep inelastic ep 
(or 7*p) scattering. The fitting process starts with an initial Qq, typically less than m^, the square of the c quark 
mass of w 2 GeV^, and individual quark and gluon initial distributions which are parameterized with pre-determincd 
shapes in x determined by a set of adjustable input parameters — given as functions of x for the chosen Qq. The 
distributions are then evolved numerically on a two-dimensional grid in x and Q^ to larger Q^ using the coupled 
integral-differential DGLAP equations 0Tal, typically in leading order (LO) and ncxt-to- leading order (NLO), and 
the results used to predict the measured quantities. The final distributions arc then determined by adjusting the 
input parameters to obtain a best fit to the data. This procedure is very indirect in the case of the gluon: the 
gluon distribution G{x, Q^) = xg{x, Q^) docs not contribute directly to the accurately determined structure function 
F2^{x,Q^), and is determined only through the quark distributions in conjunction with the evolution equations, or 
at large x, from jet data. For recent determinations of the gluon and quark distributions, see 3, 7 10). 

In the following, we will summarize our analytic method that determines the singlet structure function Fs{x, Q^) and 
G{x,Q^) directly and individually, using as input Fsq{x) = Fs{x,Qq) and Go{x) = G{x,Qq), where Qq is arbitrary, 
with the guarantee that each distribution individually satisfies the NLO coupled DGLAP equations. 

The method is extended to calculate NLO non-singlet functions, so that we can also find individual quark and 
gluon distributions analytically in terms of the starting distributions of the individual quark and gluon distributions. 



We will also give some numerical examples for non-singlet NLO valence quark distributions, comparing them to the 
MSTW [l] published NLO valence quark distributions. 

II. NLO SINGLET SECTOR 

Our approach uses an unusual application of multiple Laplace transforms [lllll2| . In this note, we use double Laplace 
transforms, first transforming the coupled DGLAP integral-differential equations into a set of coupled differential 
equations in Laplace space, and finally, into a set of coupled algebraic equations in a second Laplace space. We then 
solve the coupled algebraic equations in this second Laplace space. To obtain our final results, we must invert the 
Laplace transforms. The second transform to the algebraic equation space is analytically invertible to the space in 
which we had the coupled differential equations. The final inversionjfrom this Laplace space back to our initial space, 
must be obtained by numerical inverse Laplace transformations [l|, y] • 

We first introduce the variable v = ln(l/x) into the NLO coupled DGLAP equations. This turns them into coupled 

convolution equations in v space, which, after introducing a new variable t{Q'^,Qq) = -j- Jn^^ as{Q^)d InQ^, are 
readily Laplace transformed to obtain a set of coupled homogeneous first-order differential equations in the variable r. 
The parameters of these transformed equations are known functions of s, the Laplace-space variable. These equations 
are then Laplace transformed a second time, essentially transforming the variable r of the coupled differential 
equations into a new Laplace variable f/, with the resulting equations being coupled algebraic equations in U — with 
s again being a parameter — which are then solved iteratively. These solutions, in s and U, are analytically Laplace 
inverted back to variables s and r. Using fast and accurate numerical inverse Laplace transform algorithms [llHi ^^ 
transform the solutions back into v space, and, finally, into Bjorken x-space to obtain Fs{x, Q^) = Ts{Fso{x), Go{x)) 
and G{x,Q^) — Q{Fso{x),Go{x)), where the functions J^ and Q are determined by the splitting functions in the 
DGLAP equations. 

A similar method was used in an earlier paper jl3| in which we obtained the decoupled solutions in LO for both 
the singlet and the non-singlet sector, using only one Laplace transform. The r dependence in that case was trivial, 
and the decoupled equations could be solved directly. The extra Laplace transform that appears in the present work 
is necessitated by the nontrivial dependence of the NLO terms on r. 

Our method is readily generalized to all orders in the strong couplin g co nstant, but for brevity we limit ourselves 
to NLO in this paper. We write the coupled NLO DGLAP equations [111, [l3l schematically, using the convolution 
symbol (g), as 

'" ^^' (X, Q') = F.® (f,", + ^^P,\) {X. Q') + G® [pl + ^4^n',) (X, Q% (1) 



a,(g2)ainQ2 



a,(g2)91nQ2 



(2) 



The P°^(a;), P^g{x), Pgq{x) and Pgg{x) used in Eq. ^ and Eq. © are the LO singlet splitting function and the 
Pgg{x), Pgg{x), Pgq{x) snid P^ g [x) Av B t^B NLO siug^t splitting functions, with as{Q^) the NLO running strong 
coupling constant. It is standard procedure to construct as(Q^) assuming three massless quarks, u, d and s, below 
the c-quark threshold, adjusting the QCD parameter A at each successive threshold in Q"^ , i,e, at Q^ = M^ and M^, 
so that as{Q^) remains continuous when the number of quarks changes as heavy c and b quarks begin to contribute. 

Introducing the variable changes 

q2 

V = \n{l/x), w = ln{l/z), T{Q\Ql) = ^f a,{Q'^)d\nQ'\ (3) 

^^ JQl 

and the notation 

F,{v,t) ^ F,{e-\Q% G{v,T)^G{e-\Q^), (4) 



we rewrite the above DGLAP equations in terms of the convolution integrals 



OF 



dG , 



Fs{w,t) ( H,,{v -w) + ^K(^ -w)]dw 
+ I G{w, r) [ H,g{v -w) + '^H^iv - w) ] dw, 



— (w,t) = / Fs{w,T)[Hgq{v-w) + —^^H^^{v-w)]dw 



An 
+ 1 Giw,r)[H,,iv-w) + ^Hlg{v^w)]dw, 



where 



F,(i>,r) EE F,(e-^T), Giv,T) = Gie-^r), 



(5) 



(6) 



(7) 
(8) 
(9) 



The splitting functions Pqqi^e'"") and Pgg{e~^) in the integrals Eq. ([5|) and Eq. ([5]) involve the distribution 1/(1 — 
e~")+. The integrals involving this term can be transformed to expressions that involve derivatives of Fs or G 
without the appearance of the singular factor 1/(1 — e^^), for example, to integrals of the form /„ dFs{w, T)/dw ln[l — 
g-(«-iD)j^^ and /q" dG{w,T) / dw\n\l — e~'^^~^'']dw. After this change, all of the integrals in Eqs. (O and (|6]) are normal 
convolutions. By making a Laplace transform in u, we can factor these integrals, since the Laplace transform of a 
convolution is the product of the Laplace transform of the factors, i.e.. 



F[w\H[v — w\ dw; s 



C[F[v];s]xC[H[v];s]. 



Defining the Laplace transforms 

f{s,T) = C F,{v,t);s , g{s,T) = C[G{v,T);s] 
and noting that 



£ 



dw 



{w,t);s 



s/(5,r). 



£ 



dG 
dw 



iw,T);s 



sg(s,T), 



(10) 



[ID 



(12) 



we can factor the Laplace transforms of Eq. ^ and Eq. (|6]) into two coupled ordinary first order differential equations 
in the variable r in Laplace space s with r-dependent coefficients. These can be written as 



(.,r) = (<^fis) + 



5/ 
dr 



|(^,r) = ('5^0^) + 



".s(t), 



4^ <^r"{s)) fis,r) + [efis) + ^er°(.)) ,(.,r), (13) 

<&;--(«) 1 g{s,r) + ( e^'^(.) + ri^e--^(s) ) /(s,r), (14) 
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An 



where we recall that the the Q^ dependence is through the function r, i.e., t{Q'^,Qq) = -^ /„2 as{Q'^) d InQ'^. 
The LO coefficients ^-^^ and 6-^° are given by HI] 



'j>n^) = 4j(^+^+2(^(.+i)+,.)), 



s + 2 5 + 3/ ' 



1 



^rn, , 33-27if „ /I 2 1 



8 /2 



0."°(^) = ^ ( - - 



1 



3 \s s+1 s+2 



(15) 
(16) 
(17) 
(18) 



Here ip{x) in Eqs. ^TE\\ and (|17p is the digamma function and je = 0.5772156 ... is Eulcr's constant, quantities that 
are introduced in the Laplace transform of the LO terms involving the distribution 1/(1 — e~")+ discussed above. 
The evaluation of the NLO coefficients is straightforward, but too lengthy to be included in this note, and will be 
given, in the future, when we make numerical evaluations of Fs{x, Q^) and G{x, Q^) in NLO. 

In the case of LO, the r dependence of the equations is trivial, and the equations can be solved simply [13|, as 
already noted. The extra explicit dependence of the NLO terms on the right-hands of these equations on r prevents 
a similar construction here. In order to decouple and solve Eq. (J13p and Eq. (J14p , we Laplace transform them a second 
time — this time with respect to the variable r — into U space, i.e., we let 



T{s,U) = C[fis,T);U], 



C 



5/ 
dr 



(s,t);[/ 



= C/^(,s,C/)-/o(s), 



g{s,U) = C[gis,T);U], 
dg 



d- 



-{s,t)-U 



= Ug{s,U)-go{s), 



(19) 
(20) 



where now s is simply a parameter in U space. 

In U space, we now write the final desired coupled algebraic equations for J^(s, L'^) and g{s, U) as 



C/J-(s,f/)-/o(s) 



Ugis,U) ~ gois) 



$ 



LO/ 



asir) 



(s) J-(., U) + ^^^^\s)C[^^fis, r); U] 



An 



+e^°(s)a(s,c/) 



e 



NLO 
f 



{s)C[ 



in 



9is,r);Ul 



$iC>( 



s)g{s,u) 



$f^°(.)/:[^3(.,r);C/] 
r"s(T) 



-0^O(.)^(s, U) + ef^^{s)C['-^fis, r); [/]. 



(21) 



(22) 



For brevity, we replace the NLO "iJ by a{T) in Ec^. (PT|) and Eq. 



. We can numerically show that an excellent 



approximation to a{T) = — ^^, accurate to a few parts in 10'', is given by the expression 



a{r) 



flQ + flie" 



-biT 



(23) 



where the constants qq, oi, 6i are found by a least squares fit to a{T). We note in passing that this approximation is 
inspired by the fact that in LO, ag^LoiT) is exactly given by the form as^LoiQo)^''^ ■ 

Using the value of a^r) given by Eg. pB]). we can write the Laplace transforms C[ "'^J' f{s,T); U] 
£[^f^g(s, r); U] needed in Eq. ^ and Eq. ^ as 



and 



>c[^/(-,-);f/] 



^ajJ'(s,C/ + 6j), 






1 

Y^a,g{s,U + hj), 



(24) 



where &o = is understood in Eq. (^5)) and Eq. ([M)) . 
After introducing the simplifying notation 



$/(s)^$^-(s) 



ef{s) 



we can finally rewrite Eq. (|2T|) and Eq. ([22|) as 



ao$f^°(.), 
«oef^O(.), 



[U-<i>f{s)]F{s,U)^Qf{s)g{s,U) 

-e<,(.s)^(s, u) + [u- $3(.s)] g{s, u) 



foi-s) 
9o{s) 



- ^LOt 
- 9 
LO/ 



ao$^^° 



%{s)^<^f{s) 



63(5) = e^^(,s) + floe 



ai [$f^°(s)^(s,t/ + fei)- 
ai[ef^O(s)^(s,C/ + 6i) 



(sj, 


(25j 


^°(^), 


(26) 


0f^O(s)g(s,[/ + 6i)], 


(27) 


<^°(s)e(s,C/ + 6i)], 


(28) 



which we will solve iteratively, using ai as an expansion parameter. 

We note that the <I>'s and 0's, as defined above, contain both LO and NLO terms. We further point out that 
Eq. (|27p and Eq. (|28p are completely symmetric under the simultaneous transformations / ^ g and J- -^^ g. We 
finally remark that ai, the NLO expansion parameter in our iterative solution of Eqs. \2b\ and ([26]), is quite small: 



ai = 0.025, 61 = 10.7 for M^ < Q^ < M^ GeV^ and oi = 0.017, h = 8.63 for M^ < Q^ < 10^ GeV^, with the ao 
terms in Eq. (j25p and Eq. (|26l) being positive and about an order of magnitude smaller than the oi terms. 



We next consider the simple solutions to Eq. (l27l) and Eq. (|28p . called ^-"1(5, U) and Gi{s, U), that result from setting 
Oi = 0, i.e., the equations 



whose solutions are 



[u~'S>f{s)]Ms,u)-efis)g,iu) = Ms), 

-Q^{s)J^,(U) + [U-^gis)]g,{s,U) = gois), 



Fi{s,U) = [U-<!>g{s)]fois)/D{U,s)+efis)go{s)/D{U,s), 
Qi{s,U) = [U-'!>f{s)]go{s)/D{U,s) + eg{s)Ms)/D{U,s). 



(29) 
(30) 



(31) 
(32) 



The denominator D{U, s) in Eqs. ([31]) and ([32]) is just the determinant of the coefficients of J^{s, U) and Q{s, U) in 
Eqs. dUl) and 



(33) 



D{U,s) ^ $/(s)$3(,s)- 6/(5)63(5) -[$/(s) + $,(s)][/ + [/2 

^ {u-\ (*/(s) + '^ai^)) - \RiA {u-\ (*/(«) + *9(s)) + \r{s) 



where R{s) = yi^fis) — ^gis)) + 46/(s)6g(s). The zeros of D{U,s) lead to simple poles in Ti and Qi in the U 
plane. These functions have no other singularities, and decrease as 1/|C^| for \U\ — >■ 00. The inverse Laplace transforms 
of J-'iis, U) and Qi{s, U), denoted by /i(s, r) and gi{s, r), are therefore well defined and simple to calculate. We will 
write them as 



/i(s,t) = kff,{s,T)fa{s) + kfg,{s,T)go{s), gi{s,T) = kgg,{s,T)go{s) + kgf,{s,T)fn{s), 
where the coefficient functions in the solution are 

cosh(Ii?(.))+!^^^i^||M(^^ 



^x(^,(,H^„(,)) 2sinh(§fl(g)) 



kfgi{s,T) 



Mr) 



kgf,{s,T) = e*(*/('*)+*«('^)) 



Ris) 
cosh (— i?(s) 



Ofis), 

sinh {^R{s)) 



($/(s)-$^(s)) 



2 sinh (§i?(s)) 
W) 



63(5) 



(34) 

(35) 
(36) 

(37) 
(38) 



We comment that this solution has small NL6 terms in it, arising from the ag term in Eq. (j25p and Eq. (j26p . However, 
it is identical in form to the L6 solution we gave in Ref. |13j . and reduces to it if we set ao = 0. 

We next construct an iterative solution to Eqs. ([?7)) and (P5|) for T and Q. We start by substituting the known 
functions J-i and Gi for J-' and G on the right hand sides of the equations, and then re-solve the equations to obtain the 
next approximations J-2 and G2 for J^ and G, and then repeat the process. For the first step, we make the replacements 

J^is,U + bi)^Ti{s,U + h), Gis,U + bi)^Giis,U + bi), (39) 

on the right-hand sides of Eqs. p7|) and (p8|) to obtain our first iterative equations for J-2(s, U) and ^2(5, [/), 

[(7 - $/(.s)] J-2(s, (7) - 6/(5)^2(5, [/) = /o(s) + ai [$f ^°(.s)J-i(.s, C/ + 61) + 6^^° (5)^1(5, C/ + 61)] , (40) 
-6g(s)J-2(s,[/) + [C/-$g(s)]a2(s,f/) = 5o(s)+ai[6f^O(s)J-i(s,C/ + fei) + $f^O(s)^i(s,f/ + 6i)]. (41) 

The functions Fi{s,U) and ^i(s,C/) are given analytically by Eg. (pij) and Eg. p2p . respectively, so that we know 
them at the argument (s, U + 61), needed in the right hand sides of our iterative equations. 

Since the functions on the right-hand sides of Eqs. PO)) and (|4ip are known, and the left-hand sides have the same 
structure as Eqs. ((29| and ((30)) . their solutions can be obtained by the substitutions 

(42) 
(43) 



Us) ^ /o(.s) + ai[$f^O(s)^i(s,;7 + 6i) + 6f^O(s)gi(s,t/ + 6i)], 
5o(s) ^ go(s)-Kai[6f^°(s)J-i(s,;7 + 6i) + $f^O(s)ai(,s,[/ + &i)], 



on the right hand sides of Eqs. (|3T1) and (|32p . The leading fo{s) and 5o(s) reproduce Ti{s, U) and ^i(s, C/). The added 
terms, proportional to the expansion parameter oi, are more complicated expressions that are rational functions in [/, 
whose numerators are a second-order polynomial and whose denominators are the factorable product D{U, s)D{U + 
bi, s). The functions J-2 and Q2 therefore have an extra pole in U that is displaced along the real axis from the poles of 
J-"i and Qi by the amount bi. Since this is the only new singularity, whose terms decrease at least as rapidly as l/[/^ 
for [/ — !► 00, the overall behavior of the iterated solution decreases at least as rapidly as 1/|C/|, so that the inverse 
Laplace transforms needed can be calculated analytically. 

We can again write the results for the inverse transforms /2(s, r) and 52(s, t) in terms of the initial distributions 
/o(s) and go{s) as in Eq. ((34)) . but with the coefficient functions k now sums of the original expressions in Eqs. ((35)) - 
((55)) and terms that depend linearly on the coefficient Oi in Eq. ((^5)) as well as on s and r. The coefficient og has been 
incorporated into the definitions of the $'s and 6's in Eqs. ()25p and ((^5)). so it does not appear explicitly. 

Continuing, we find the k^^ iterated solution to Eqs. ((?7)) and ()28p for J-'{s, U) and Q{s, U) by making the substi- 
tutions 

/o(.s) ^ /o(.s) + ai [$f ^O(s) J-fc(s, [/ + &i) + ef ^0(s)gfc(s, U + 61)] , (44) 

5o(s) ^ /o(s) + ai [<i>f ^°(s)efe(s, [/ + 61) + ef ^0(s)J-fe(s, [/ + 61)] (45) 

in the right-hand sides of Eqs. ((3T)) and ((32p and replacing J^i(s, C/) and 5i(s, t/) on the left-hand side by Fk+i{s, U) 
and Gk+i{s, U). The resulting expressions for J-fc+i and ^fc+i add new terms proportional to aj^, which again are ratio- 
nal functions of U, with denominator of higher power in U than the numerator. All the terms in Tk+i{s, U) , Qk+i (s, U) 
decrease at least as rapidly as 1/|[/| for \U\ — >■ 00, and the only singularities are poles at known locations, so we can 
again calculate the Laplace inversion from U space to r space analytically. At each stage, we can write the inverse 
transforms as 

f{s,T) = kff{ai,bi,s,T)fo{s) + kfg{ai,bi,s,T)go{s), (46) 

g{s,T) = kgg{ai,bi,s,T)go{s) + kgf{ai,bi,s,T)fo{s), (47) 

with the functions fc(ai, 61, s, r) expressed as power series in the NLO expansion parameter oi whose coefficients are 
analytic functions of s and r. These expressions rapidly become too complicated and too lengthy to reproduce here, 
but are easily calculated using a program such as Mathematica |14) |. 

After numerical Laplace inversion [ij, [2| of the fc's from s to w space, suppressing their explicit dependence on ai 
and 61, we define their Laplace inverses as 

Kff{v,t) = £"-^[fc//(s, t);w], Kfg{v.,t) = C~^[kfg{s,T);v\, (48) 

Kgg{v.t) = C-^[kgg{s,T);v], KGF{v.T)=C-^[kgf{s,T)-v], (49) 

SO that wc can write the decoupled solutions in (i?, r) space as the convolutions 

pV pV 

F,{v,Q^) = / KFF{v~w,T{Q\Ql))Fso{w)dw+ Kfg{v - w,T{Q'\Ql))Go{w)dw, (50) 

Jo JQ 

2\ — / j^ /„. „,, ^tn!i /n2NN A_ ,'„,,\ J ,, I / 7>' /„. „,, ^(rt'2 /o2 



G{v,Q') = / KGG{v-w,T{Q\Qi))Go{w)dw+ Kgf{v ~ w,T{Q\Qf,))F,o{w)dw. (51) 

Jo Jo 

Finally, recalling that v = ln(l/a;), we can transform the above solutions back into the usual space, Bjorken-x 
and virtuality Q^, enabling us to write the NLO decoupled solutions, Fs{x,Q'^) and G{x,Q'^), which require only a 
knowledge of Fso{x) and G{x) at Qq, where evolution is started. 

In order to insure continuity across the boundaries Q^ — M^ and M^, we first evolve from Qq (where, e.g., Qq = 1 
GeV^ for the MSTW group [3| ) to M^ and use our evolved values of Fgo (w) and Go (v) for a new starting values of 
Fso{v) and Go{v). We then evolve to M^, repeating the process, thus insuring continuity of Fs{x, Q^) and Gs{x, Q^) 
at the boundaries where nf changes. 

III. NON-SINGLET SECTOR 

For non-singlet distributions Fns{x, Q^), such as for valence quarks, -Dvai = x (d{x, Q^) — d{x, Q^)) — the difference 
between quark distributions — we can schematically write the logarithmic derivative of F„s as the convolution of 



Fns{x,Q'^) with the non-singlet sphtting functions, PniP'^^ix) and Pj^^'^'"*(a;), for LO and NLO, respectively, (using 
the convolution symbol Cg)), i.e.. 



47r dEns ,__ ^2 



(.T,g2) = F„ 



pLO,ns , (^s(T) pNLO.ns 
11 47r 99 



a,(g2)ain(Q2) 
After changing to the variable v = hi(l/a;) and the variable r, we write 



(x,g2). (52) 



or ./n 



An 



P^O,ns^^_^^^<^pNLO.ns^^,_^^ 



dw. (53) 



The comments that we made in Sec. II about integrals that involve the distribution 1/(1 — e^^)+ also apply here. 

Going to Laplace space s, we obtain a linear differential equation in r for the transform /„s(s,r). This has the 
simple solution 



/„,(s,r) = e-*'"(-^)/„,o(5), $„.(.) ^ $,i°(.) + ^$^,^°(s), (54) 



where 



and 



47r Jo (47r)-^ Jq2 



(55) 



We note that in LO, $ns(s) = <I'^'^(s), where $^^^(3) has been written out explicitly in Eq. (fT5|) . Again, the evaluation 
of ^ns^'^i^) i^ straightforward, but too lengthy to be shown here. 

We can find any non-singlet solution, Fns{x, Q^), by using the non-singlet kernel Kns{v) = £^^ [e'^*"^^'*^; u] in the 
Laplace convolution relation 

Fns{v,T)^ / Kns{v -W,T)Fnso{w)dw, (57) 

Jo 

and then returning to (a;, Q"^) space. 

In order to insure the continuity of Fns{x, Q^) where Uf changes, we renormalize the starting values Fnso{v) at the 
boundaries M^ and M^, as described previously in a similar context for singlet distributions. 

A. Comparison of non-singlet theory with NLO MSTW non-singlet valence quark distributions 

As an example of the application of this technique, we will compare two x-space non-singlet valence quark dis- 
tribution functions i^„s(a;, Q^) calculated from Eq. (|57|) with the published MSTW values Q- In Eq. ()57p . we use 
Ql = I GeV^ the MSTW starting value for evolution, to construct Fn^oiv) from the pubhshed NLO MSTW Q 
quark distributions. We use the MSTW values Mc = 1.40 GeV, Mb = 4.75 GeV, together with the MSTW NLO 
definition of as(g2), adjusted to be continuous at the boundaries Q^ = M^ and M^, with as(l GeV^) = 0.49128 and 
a^(M|) :=: 0.12018 g. 

1. The NLO non-singlet d quark valence distribution Dvai ~ x (d{x,Q^) — d{x,Q^)) 

In Fig. [1] we show the results obtained by evolving the non-singlet d quark valence distribution, Dvai = 
x{d{x,Q'^)-d{x,Q'^)), from Ql = 1 GeV^, for Q^ = 5 ^20, 100 and M^ GeV^. The published MSTW |3 curves 
are: Q^ = 5 QeV^, solid blue; Q^ = 20 GeV^, dashed green; Q^ = iqo GeV^, dot dashed red;_Q2 = m^ GeV^, large 
dashed black. The dots are our evolution results for NLO non-singlet Z?vai = x {d{x,Q^) — d{x,Q^)) from Eq. ((57|) 
(converted to x-space), using the NLO MSTW values for Fnso{x), where Qq = 1 GeV^. Let us define the fractional 
error, frac. err. = 1 — Dvai (calculated) /I?vai (MSTW), at a; = 0.135, a point near the peaks of the curves in Fig. [T] 
The reproduction of the published MSTW data is excellent. We find that at Q^ = 5 GeV^, frac. err. = -0.004 and 
at Q2 = M^, frac. err. = -h0.004. 
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FIG. 1: The NLO MSTW 0] non-singiet valence distribution, D^^i = x {d{x,Q'^) - d{x,Q'^)), for Q^ ^ 5 ,20, 100 and M^ 
GeV^ The pubhshed MSTW H curves are: g^ = 5 GeV^ sohd blue; Q^ = 20 GeV^ dashed green; Q^ = 100 GeV^ dot dashed 
red; Q^ = Mf GeV^, large dashed black. The dots are the evolution results for NLO non-singlet Dvai ~ x {d{x, Q^) — d{x, Q^)) 
from Eq. ([57)l (converted to a;-space), using the NLO MSTW values for Fnso{x), where Qo = l GeV^. 




FIG. 2: The NLO MSTW [3 non-singlet valence distribution, t/vai = x {u{x,Q^) - u{x,Q^)), for Q^ = 5 ,20, 100 and Mf 
GeV^ The published MSTW [3 curves are: Q^ = 5 GeV^ solid blue; Q^ = 20 GeV^ dashed green; Q^ = 100 GeV^ dot 
dashed red; Q^ = Aff GeV^, large dashed black. The dots are the evolution results for NLO non-singlet valence distribution, 
f^vai = X (u{x,Q^) — u{x,Q^)^ from Eq. (|57|) (converted to s-space), using the NLO MSTW values for F„so{x), where Qo = l 
GeV^ 

2. The NLO non-smglet u quark valence distribution f/vai = x (u{x,Q^) — u{x,Q^)) 



In Fig. [21 we show the results obtained by evolving the non-singlet distribution valence distribution for the u quark, 
Uy^x, g2) = X {u{x, Q2) _ u(a;^ g2)^ ^^.^^ q2 ^ ^ Q^y2^ fQ^ Q2 = 5 , 20, 100 and M^ GeV^. The published MSTW 

3 curves are: Q^ ^ 5 GeV^, solid blue; Q'^ = 20 GcV^, dashed green; Q^ = 100 GcV^, dot dashed red; Q^ = M^ 
GeV^, large dashed black. The dots are our evolution results for NLO non-singlet u quark valence distribution 
Uvuiix, Q^) = X (u{x, Q"^) — u{x, Q^)) from Eq. ([57)1 (converted to x-space), using the NLO MSTW values for F„so(a;), 



where Qq — 1 GcV^. Again, the reproduction of the published MSTW data is excellent. The fractional errors at 
X = 0.135, defined in Section [mXTl are: frac. err. = -0.003 at Q^ = 5 GeV^ and frac. err. = +0.004 at Q^ = m^ 
GeV2. 

IV. CONCLUSIONS 

For the singlet sector of pQCD, we have solved the coupled NLO DGLAP equations and found NLO decoupled 
analytic solutions for Fs{x,Q^) and G{x,Q^), an extension of our earlier work for LO [ij]. All that is required is 
knowledge of the initial distributions Fsq{x) and Go(a;), at Q^ — Qo, where Qq is the starting value for the evolution. 
For the non-singlet sector, we have successfully solved the NLO evolution equation for Fns{x,Q'^), again in terms of 
Fnso{x), the value oi Fns{x,Q^) at Qq. We illustrated this numerically for NLO, calculating the non-singlet valence 
quark distributions C/vai = x (u{x, Q^) — u{x, Q^)) and Z?vai = x [d{x, Q^) — d{x, Q^)) for a very large range of x and 
Q^, in excellent agreement with the NLO published MSTW Q values. We note that these techniques can be extended 
to arbitrary order in the strong coupling constant as{Q^), for both the singlet and non-singlet sector. 

The results presented here are basically analytic, thus eliminating the need for simultaneous numerical solutions of 
the singlet and non-singlet DGLAP equations on a two-dimensional lattice in x and Q^. They provide new tools for 
studying pQCD; for example, they can be used to examine directly the sensitivity of an individual evolved distribution 
to the assumed shapes of its starting distribution. In the future, we hope to apply these techniques to a global fit of 
experimental F2^{x,Q'^) data to determine in LO, a gluon starting distribution, and in NLO, approximate Fs and 
gluon starting distributions. Since these starting distributions will be determined by experimental data, they will be 
free of predetermined shape hypotheses; this will allow us to find new gluon distributions that are critically needed 
for the interpretation of results from the Large Hadron Collider. 
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